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Abstract 

We propose a scheme for implementing cross Kerr nonlinearity between two superconducting 
transmission line resonators (TLR) via their interaction with a coupler which is constructed by two 
superconducting charge qubits connected to each other via a superconducting quantum interference 
device. When suitably driven, the coupler can induce very strong cross phase modulation (XPM) 
between the two TLRs due to its N-type level structure and the consequent electromagnetically 
induced transparency in its lowest states. The flexibility of our design can lead to various inter- 
TLR coupling configurations. The obtained cross Kerr coefficient is large enough to allow many 
important quantum operations in which only few photons are involved. We further show that this 
scheme is very robust against the fluctuations in solid state circuits. Our numerical calculations 
imply that the absorption and dispersion resulted from the decoherence of the coupler are very 
small compared with the strength of the proposed XPM. 
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I. INTRODUCTION 



The circuit quantum electrodynamics (QED) which employs the superconducting 

coplanar transmission line resonator (TLR) to substitute the standing-wave optical cavity 
and superconducting qubits to replace the atoms is an on-chip realization of cavity 
QED [6]. Compared with conventional optical implementations, this solid-state architecture 
offers unprecedented tunability and scalability which are leading to flexible quantum optics 
in electronic circuits. Since the strong coupling between the TLRs and the superconducting 
qubits with vacuum Rabi frequency up to three orders larger than the qubit decoherence 



and cavity decay has already been achieved Q], many important quantum information 
processes (QIP), including coupling qubits using the TLR as a data bus llO,LLl| and preparing 



the TLR Fock states 19, 12 



19. 



13] have been demonstrated. 



Stimulated by the advances on the single TLR level, recently the ideas of photon ma- 
nipulations between TLRs have been developed 15-[l7] . The motivation is to facilitate the 
future realization of scalable quantum computation. These schemes mainly consider a model 
of two TLRs coupled to an assistant tunable coupler. One can tune on and off the individ- 
ual energy transfer between the TLRs and the coupler through frequency selection and thus 



manipulate the two mode photon states. A very latest experiment have rea 



1* 



ized the NOON 



state preparation in two TLRs connected to an entanglement generator 

In optical systems, besides the linear tight-binding photon transfer, the nonlinear Kerr 
interaction between cavity modes has also been studied extensively. The cross Kerr nonlin- 
earity, or the so-called cross phase modulation (XPM), has found wide application in QIPs 



including the construction of nontrivial quantum gates 
photon states 



21 



24 



- |23| , the preparation of entangled 



341 ]. and quantum non-demolition measurement [25]. Enhancement of 



dissipation-free photon-photon interactions at the few-photon level is a fundamental chal- 
lenge in quantum optics. Since photons can hardly interact with each other, the XPM is 
often obtained by coupling two photon modes to an atomic nonlinear media. To minimize 
the dispersion and absorption, schemes of exploiting destructive quantum coherence in N 



type atoms have been proposed 



261 ] and realized 



27 



291 ] . The reported experiments are 



performed in the semiclassical region due to the very small interaction strength between the 
laser fields and 87 Rb atoms, which is often on the same level of the atomic decoherence rates. 
In order to obtain significant phase shift, the probe and control pulses contain large numbers 
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of photons. The weak coupling places a major hindrance to the further application of the 
Kerr effect on the single photon level. 

In this paper, we aim at the realization of strong cross Kerr coupling between two TLR 
modes. This work is inspired by the recent self phase modulation scheme in circuit QED sys- 
tem (41). We design a superconducting circuit which exhibits complete analogy to the N level 
XPM schemes in atomic systems. The proposed four level artificial molecule is constructed 
by two Josephson charge qubits coupled by a superconducting quantum interference device 
(SQUID). When capacitively coupled to the two TLRs, the molecule induces cross Kerr in- 
teraction between microwave photons in the two TLRs. The strong TLR-molecule coupling 
can boost the XPM strength up to several MHz. Moreover, the dispersion and absorption 
resulted from the decoherence of the coupler are estimated to be negligible compared with 
the obtained XPM strength. Our system is flexible enough to allow various inter-TLR cou- 
pling configurations. Since long lifetimes for both the TLRs and the molecule have already 
been achieved, many QIPs between the TLRs in which there are only few photons involved 
can be performed with very high fidelities. 

The paper is organized in the following manner. We first briefly describe the general 
N-type scheme of XPM in Sec. [XT] and then study its realization in circuit QED system in 
Sec. IIII1 In Sec. IIV} the influence of the decoherence on the XPM scheme is investigated in 
detail. The application of our XPM scheme and the related discussion are presented in Sec. 
IVl while the conclusion is given in Sec. I VI I 

II. XPM INDUCED BY THE N TYPE ATOM 

We start from the theoretical model of two cavity modes, labeled by their annihilation 
operators a\ and coupled to an atom with N type level structure, as shown in Fig. [TJ 
The first mode couples the |1) -H- |3) transition with strength g\ and detuning 5, while the 
second mode couples the |2) <-» |4) with strength g 2 and detuning A. In addition, a classical 
laser field is applied to drive the |2) -h- |3) transition with strength Q c . The frequency of 
the classical field is set to be in Raman resonance with the first mode. In the interaction 
picture, the Hamiltonian of the system can be written as 

H S ys H atom Hint-, (1) 
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FIG. 1: (Color online). The N scheme for cross Kerr interaction. Two cavity modes a\ and a2 
induce the transitions |1) — > |3) and |2) — > |4) of the N type atom with coupling strength gi and g 2 
respectively. The 5 and A are the corresponding detunings, while the factors 71, 72, and 73 label 
the decay rates from the upper levels to the lower levels. In addtion, an external pumping pulse 
couples 1 2) and |3) with Rabi frequency $7 C . The classical pumping is set to be in dark resonance 
with the first cavity mode so that the states |1), |2) and |3) exhibit electromagnetically induced 
transparency. 



with 



H, 



atom 



Sa 33 + A(T 44 , 



hit 



i[gi{a[vi3 ~ 0- 3 lOl) + 02(4^24 - 0-420-2) 



(2) 
(3) 



where Ojk = \j) (k\ are the atomic raising and lowering operators for j 7^ k and population 



operators for j = k. In the limit 



20-122 1 



9i 



< 1, \g 2 \ < |A| 



(4) 



the atom evolves rapidly on the time scales relevant for the cavities. After the adiabatical 
elimination of the atomic degrees of freedom, an effective Kerr photon-photon interaction 



TT ~ 

-"SJ/S — 



H, 



92 9i t t 
eff - --^-^2 0, 1 a 1 a 2 a 2 , 



(5) 



can be obtained 



26] 



The Kerr effect induced by the N type atoms has been investigated in the past few years 



19 



22 



26 



-|29|. The physics of H e ff can be interpreted in an intuitive way: Suppose initially 
the atom is in its ground state |1). In the first step we set g 2 = so that |1), |2) and |3) 
decouple from |4). When the classical pumping is added, the states |1) and |2) form a dark 
state, in which the destructive quantum interference between the |1) -H- |3) and |2) -H- |3) 



transition branches cancels the dispersion and absorption 35j of the a\ mode. The stationary 



population on the level |2) is thus determined by the photon number a\a\ as 

(a 22 ) - gfal^/ni (6) 

When the second mode is added back to couple the |3) -H- |4) transition dispersively, the 
resulting AC Stark shift has the form — g 2 a\a 2 a 22 /A = — g\g\a} x a\a} 2 a 2 l (Af^). 

Mathematically, to derive H e ff, we first write the evolution equations of the operators a\ 
and a 2 as 

= ~ Kiai ~ V2tTia lin (t), (7) 

= ^20-24 - K 2 a 2 - \[2K 2 a 2in {t), (8) 
where Ki j2 are the decay rates of the cavities, and a lt2in (t) are the input noise operators 
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331 ] . Adiabatic elimination of the atom allows us to express the atomic operators in 
terms of the mode operators. Formally this is accomplished by setting the time derivatives 
of all the <jjk to be zero. The stationary values (<Jjk)s °f the atomic operators can be 
expanded as 

oo 

(^>s = 5>^ ( 9 ) 

n=0 

where a™ k contains the multiplication of n creation or annihilation operators of the cavi- 
ties. We further use cru = 1 as the starting value and iteratively determine the remaining 
coefficients of the expansion. We thus get 

2 

(0-13)5 = -f^ a l a ^ ( 10 ) 
2 

/ \ 9i92 t 

(024>S = ~72^2 a l aia2 - C 11 ) 

Replacing CX13 and a 2 4 in Eqs. ([7]) and ([8]) by their stationary values in Eqs. f fTUl) and f fTT]) 
we obtain 

dcLi 9^9? t 1 / \ / \ 

—jr = -TT-FTJ a 2 a i a 2 - - VS/tiaa^^J, (12) 



~1T = _ ^T7^" a i a i a 2 - ^2^2 - V^2a 2in (t)- (13) 
at i/\\¥ c 

The first terms in Eq. f|12j) and ([TBI yield an effective Kerr interaction 

H e ff = — rja^ala^, (14) 

where rj = g\g\j (Af^) is the Kerr coefficient. 

The first order terms in the expansion of (013)5 anc ^ ( a ^)s which represent the linear 
dispersion and absorption of the two cavities, are vanishing. The suppression of linear 



susceptibilities is a result of the quantum phase coherence in the atom [35] . When the 
decoherence of the atom is taken into account, the phase coherence is broken and the linear 
susceptibilities become nonzero. Also, the third order susceptibility 77 is modified from real 
to complex. The derivation of Eqs. ffTUl) and pip and the discussions about the influence of 
the atomic decoherence are provided in the forthcoming sections and the Appendix. 



III. THE N SCHEME IN CIRCUIT QED 

TLR A 















Coupled CPB 

















TLR B 



FIG. 2: (Color online). Schematic circuit for XPM between two TLRs. The TLR A (blue) and 
the TLR B (red) are capacitively connected to a superconducting molecule (gray) controlled by an 
external circuit (black). 

The N scheme described in Sec. [Ill can be realized in superconducting quantum circuits. 
We consider a system of two TLRs capacitively coupled to an artificial superconducting 
molecule, as shown in Fig. |2j In the following subsections we explicitly show that all the 
required ingredients of the N level XPM proposal can be established in the circuit QED 
system. 
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A. The superconducting molecule: an artificial N type system 



The artificial molecule in Fig. [2] is constructed by two superconducting Cooper pair boxes 



(CPB) |36|, |42| coupled by a SQUID 37], as sketched in Fig. EJ Each CPB consists of a 



small superconducting island connected to the ground electrode by a symmetric SQUID 
with capacitance Cji and tunable Josephson coupling energy Eji for i = 1,2. The gate 
voltages V g i bias the corresponding qubits via the gate capacitors C g i. Finally, the CPBs are 
connected to each other by a coupling SQUID with capacitance C m and Josephson energy 
Ej m . The Hamiltonian of the molecule reads 



H = 4£ m (ni - n gl )(n 2 - n g2 ) - E Jm cos($i - $ 2 ) 
+ [E ci (ni - n gi ) 2 - E M cos $J , 



(15) 



t=l,2 



where n s i )2 = C g i j2 V^i j2 /2e denote the gate- induced charge numbers on the CPBs, $i 2 are 
the canonical conjugate variables to n 12 , E cl2 = 2e 2 Cs 2i i/(CsiCs 2 — C^J are the effective 
Cooper-pair charging energies (Ca = C g i + Cji + C m is the sum of all capacitances around 
the zth qubit), and E m = e 2 C m / '(C^iC^ — Cm) * s the capacitive coupling strength between 
the CPBs. Near the co-degeneracy point n gl = n g2 = 1/2, we can use the two- level language 



71; 



;i + a x . i )/2,cos$ i = -a zi /2, 



(16) 



for i = 1,2, to describe the molecule as 



E m [a xl a x2 - b (a zl (r z2 + a yl a y2 )} 

+ - ^ [Ebi&xi + Eji<T z i] , 



(17) 



i=l,2 



where E^ 2 = 2 [E cl)2 (1 — 2n g i j2 ) + £? m (1 — 2n s2j i)] are the effective charge biases, and b 
Ej m /4:E m . 



Being a Coulomb b 



degrees of freedom 



43 



ockade device, the coupled-CPBs is very sensitive to noise from charge 



52j. By operating the molecule at its optimal point, chosen such that 



the linear longitudinal qubits-noise coupling vanishes, we can prolong its dephasing times 
by several orders 42] . For this reason, here we concentrate on the behavior of the coupled 
CPBs at the co-degeneracy point n g \ = n g2 = 1/2, which can be verified as the optimal 
point later. Without loss of generality, we further assume that the two CPBs are identical, 
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i. e. Eji = Ej2 = Ej, Cjx = Cj2 = Cj, and Csi = Cs2 = Cs- In this situation, the 
eigenstates and eigenvalues of H can be written as 



|1) = -sin#|00) + cos0|ll), 
|2) = (-|01) + |10))/ v / 2, 
|3) = (|01) + |10))/V2, 
1 4) = cos#|00) + sin6»|ll), 

E\ = —E mN — E m b , 
E 2 = -E m (l-2b ), 

E3 = E m , 
E4 = E m N — E m bo. 

where E mN = ^ E] + E*Jl - 6 ) 2 , E = V ' E j + Eg, and 9 
aicsin(Ejb /2E mN )]/2. 

V* 



n 2 ,<b 2 




F C 



F C 



(18) 



(19) 



[arcsin (i? m /i?) + 



FIG. 3: (Color online). Schematic plot of the coupled CPBs. The crosses label the Josephson 
junctions, while the round circle in the SQUID loops represent the penetrating flux. 

The levels {|1), |2), |3), |4)} in Eq. f|T8|) have the required N configuration in Sec. |TT| 
and can be mapped to the four states in Fig. [1] one by one. In experiments, Ej is often 
larger than E m , while the maximum of bo is usually on the order of unity. We thus can 
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modify the rough energy scale of the molecule by tuning Ej. Moreover, bo can be used to 
control the details of the energy structure. When bo = 0, the N level is symmetric, i. e. 
E4 — E 2 = E3 — Ei. Tuning 6 to non-zero value can break the symmetry. The influences 
of b on the level spacings E i2 = Ei — E 2 , E 3 i = E 3 — Ei, and E 32 = E 3 — E 2 are shown 
in Fig. HJ In this calculation, we set the parameters of the coupled CPBs as Ej/2tt = 20 
GHz and E m /2ir = 5 GHz {39), As plotted in Fig. H] the level spacings go through very 
large range when bo varies in the region 60 £ [0,0.8]. Therefore, with the tunability of 60 
and Ej, the level splits of the molecule can be modulated at will. We further notice that the 
co-degeneracy bias point n g i = n g2 = 1/2 remains to be the optimal point for the molecule 
during the tuning of b . The operators a x i and a x2 , through which the system couples to 
the charge noise, have the form 



&xl 






— sin^ 


cos 4 





— sinc/i 








cos 


COS (f) 








sin 





cos 6 


sin (b 






(20) 



cr X 2 



(21) 



sin0 cos0 
sin — cos q 
cos (ft sin0 
— cos (f> sin (p 

where <fi = 9 + ir/4. We thus verify that the diagonal entries of a x i and a x2 are all zero, 
which indicates that the linear longitudinal dephasing vanishes and the molecule is subject 
only to the second order, quadratic dephasing. 

Our design can be regarded as an innovation of the N type system proposed in Ref. 



4=11 ] . where Rebic et. al have proposed to build an N level molecule by two capacitively 
coupled CPBs. In that situation, the N level structure is always symmetric on the optimal 
point. To get asymmetry, additionally DC charge bias was needed and only moderate 
I £42 — -E31I can be obtained. Moreover, since the main decoherence source of the charge 
based superconducting circuits is the low-frequency charge noise, the DC bias in previous 
SPM scheme results longitudinal dephasing which severely damages the phase coherence of 
the molecule. Therefore, compared with previous design, our introduction of the coupling 
SQUID can offer more tunability and robustness to the artificial molecule. 
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FIG. 4: (Color online). The level spacings versus bo. The solid, dashed, and round lines represent 
E31, E42, and I?32, respectively. 

B. The classical pulse 

We proceed to show how to implement the classical pumping connecting |2) and |3). 
Penetrating the external flux into the SQUID loops of the CPBs can couple microwave 
pulses with the molecule through the operators o z \ and a z2 , which have the matrix form 



O'zl 



cos2# -sin2# 

0-10 
0-10 

•sin 26 cos 26 



(22) 



Oz2 



cos 26 
1 
1 



sin 26 





(23) 



-sin 20 cos 20 

The matrix elements of o z \ and a z2 support the transitions inside the Hilbert subspaces 
span{|l), |4)} and span{|2), |3)}. Therefore, we can apply an AC flux pulse to the SQUID 
loop of the first qubit to modulate Ej± with amplitude Qex and frequency u> p . When uj p 
is tuned close to E32, the |2) <H- |3) transition is effectively induced while the |1) <H- |4) 
transition can be neglected due to frequency mismatch. To further suppress the unwanted 
|1) <H- |4) transition, we can synchronizely modulate Ej 2 with the same amplitude and 
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opposite phase as that of Ej%. In this way we get a pumping Hamiltonian 

Hp Um p = Qex cos uj P t(cr zl - a z2 ) (24) 
^ -tt Ex [exp (iupt) |2) (3| + exp (-iuj P t) |3) (2|] 

which establishes the pure |2) <->■ |3) transition. 

C. Coupling between the TLRs and the molecule 

We consider the capacitive coupling between the TLRs and the molecule. As shown in 
Fig. |2l The TLR A (blue) has length L A and fullwave frequency u A , while the TLR B (red) 
has length L B and fullwave frequency u B . The Hamiltonian of the two individual TLRs are 

Hcavity = UJAa'a + u B b\ (25) 

where u A , B = 2ir/ (^La^VFcJ are the eigenfrequencies, F and c being the inductance and 
capacitance per unit length, La,b being the length of the two TLRs, and a, b are the 
annihilation operators of the fullwave modes, respectively. The quantized voltages of the 
TLRs are 

V A (x) = y/u A /L A c((J + a) cos(2nx/ L A ) , x E [-L A /2, L A /2] , (26) 
V B {x) = y/u B /L B c{tf + b) cos(2ttx/L b ), x e [-L B /2, L B /2] . 

When the two TLRs are capacitively connected with the molecule, the interaction Hamil- 
tonian has the general form 

H cc = V x [CaiVa(xai) + C B1 V B (x B1 )} + V 2 [C A 2V A (x A2 ) + C B2 V B (x B2 )} , (27) 

where the C A j-, C B j for j = 1,2 are the coupling capacitance between the jth qubit and 
the TLRs, the x A j, x B j for j — 1, 2 are the locations of the coupling capacitor, and V\ y2 = 
[C-£(Jxi,2 + C m a x2j i\ j (C| — C^J are the quantized voltages of the first and second CPB. We 
further expand H cc as 

H cc = h A i(a) + a)a x i + h A2 (a) + a)a x2 (28) 
+h B1 (b t + b)a x i + h B2 (b ] + b)a x2 
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where h A ,B;i,2 are the coupling factors which can be written as 

h A ,B;l,2 = J J -7^2 77JTpA,B;l,2COs( )C E2 ,1 (29) 

V l a,bC (C e - G^J L a ,b 

+ C j4 ,S;2,l cos( — )6 m J. 

J^A,B 

The operators er^i and cr^ have matrix elements which induce the transitions between 
subspace span{|l), |4)} and span{|2), |3)}. We can thus tune the Ej and b so that E± 2 
and £"31 are close to the cavity mode frequencies ua,b- 111 this case the TLRs can effectively 
only induce the |1) <H- |3) and |2) •<-)■ |4) transitions. The remaining |1) -h- |2) and |3) -H- |4) 
transitions are suppressed due to frequency selection. With the rotating wave approximation, 
H cc finally reads 

H cc = gAi{cri3(J + o"3i a) + 9Bi(cri3^ + 031 &) (30) 
+fiU2(cT24a f + cr 42 a) + ^2(^24^ + cr 42 6), 

where the coupling factors (?a,_b ; i,2 are 

9A,B;1 = COSip (h AtB;1 + h A ,B;2) , (31) 
9A,B;2 = COS if (h A:B ;l - h A ,B;2) ■ 



D. XPM in circuit QED 

H cc in Eq. (1301 has a more general form than the atom-photon coupling terms in Eq. ([3]) . 
The flexibility of superconducting devices allow us to choose suitable values and locations 
of the coupling capacitors in order to get desired coupling configurations. The setting of the 
locations of the capacitors does not require to put the CPBs in distant places, because the 
TLRs can be fabricated in a zig-zag form. A trivial case is that Cbi = Cb2 = 0, i. e. the 
molecule is connected only to the TLR A, which results the SPM of circuit QED [41]. To 
establish the XPM described in Sec. DTI we se t the coupling capacitors to have capacitances 
Cbi = Cb2 and Cm = C A i. In addition, the location of the capacitors are selected as 
%A\ = x A 2 = and xbi = Lb/2 — xb2 = L B /8. In this way, h A1 = h A2 , h B \ = —hs2 and the 
resulting H cc is given by 

H cc = gM^a 1 + cma) + ^2(^24^ + v^b). (32) 
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Up to a trivial % factor, H cc in Eq. ( 132|) has the same atom-photon coupling form as that 
of Hi n t in Eq- ©, i. e. the TLRs A and B replace the role of the cavities a\ and a 2 , 
respectively. 

Combining Eqs. ( TTBj) . ( 12~4|) . and ( 13"i2"j) . We develop all the required elements of the N level 
XPM scheme in the circuit QED system. When we tune the classical pumping to be in dark 
resonance with the TLR A, an effective Kerr interaction hamiltonian 

Heff = - 9 -f^LaWab, (33) 

^ lL Ex 

can be induced between the two TLRs. 

The Kerr coefficient % 3 = g\ x g 2 m l (AQ Ex ) can be estimated based on the reported ex- 
periments. Since the coupling strength between a TLR and a qubit as large as 300 MHz 
has been achieved, we set gi/2ii = (72/2vr = 300 GHz. In addition, the classical pumping 
strength Qex and the detuning A = E i2 -ub can be chosen as Qex/Zk = 1-5 GHz and 
A/2vr = 1.5 GHz to fulfill the adiabatical condition in Eq. ®. We then get Xs/^ — 2-5 
MHz, which has already exceeded the observed Kerr strength obtained by exploiting a large 



Josephson junction connected with the TLR in the recent papers 



4J-|46|. 



IV. THE ROLE OF ATOMIC DECOHERENCE IN XPM, 

We discuss the influence of molecule's decoherence on the XPM. As shown in Fig. [TJ in 
the first step we consider only the decay processes |3) —> |1), |3) — > |2), and |4) — > |2), with 
decay rates 71, 72, and 73, respectively. In this case the polarizations (<Ji3) s and (0-24)3 have 
the form 

2 

{*n) s = - {l3 + tA)n f 2aia 2 , (34) 
2 

/ \ _ 9i92 f 

(73 + i A) Q 2 C 

which results a complex Kerr coefficient 

We notice that X3 does not depend on either 71 or 72. In addition, the first order terms of 
(o" 13 ) 5 and ((J24) are still zero. These two effect can be explained by the EIT in the lowest 
levels of the molecule. The destructive interference between the |3) -H- |1) and |3) o |2) 

13 



transition branches cancels the first order terms of (013)3 and (024)5 and thus the linear 
susceptibilities of the TLRs. Moreover, the coherent population trapping significantly sup- 
press the population (033)3. As a result, the decay channels |3) <-> |1) and |3) <-> |2) become 
irrelevant. Therefore, we use the relative decay rate R define by R = {^(xs) /^(Xs)\ to 
characterize the molecule- induced cavity decay. Since 73/2^ has been pushed to the order 
of 0.5 MHz in recent experiments 8|, while A can be usually set to be on the order of GHz, 
R = 73/A is estimated to be in the range [10~ 3 , 10~ 4 ]. 




Y 5 I2k (GHz) 



FIG. 5: (Color online). The Kerr nonlinearity x'3 = ^te) versus &ex and 75. 




FIG. 6: (Color online). The linear dispersion factor X1/X3 versus Qex and 75. 

We further take the decay |2) — > |1) and dephasing between |2) and |1) into account. 
The corresponding decoherence rates are 74 and 7^. These two processes can break the 
phase coherence between |2) and |1), which is essential in the previous suppression of linear 
dispersion and absorption. In this circumstance the first order of (013)5 emerges. It takes 
the form 

°13 = -1fc9Aia/B, (36) 
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FIG. 7: (Color online). The linear absorption factor X1/X3 versus Q.Ex and 75. 




FIG. 8: (Color online). The nonlinear absorption factor 9(x3)/X3 versus Qex an d 75- 

where 75 = 74 + 7^ and B = Q Ex + 7175 + 7275 + iSjs. The resulting linear susceptibility is 

Xi = g 2 A1 ls/B. (37) 

the real part x'i °f Xi which describes the decoherence induced dispersion of the TLR A is 
approximately g 2 Al l$ / & Ex . With 7 5 /2vr = 0.5 MHz, x x /2>k is lower than 0.01 MHz. The 
imaginary part x'l °f Xii which describes the decoherence induced dissipation of the TLR 
A, can be approximated as x'{ — QaiI^Is/^ex- ^ i s even smaller than x'i by several orders. 
The analytic form of X3 m the presence of the decoherence |2) — > |1) and the dephasing 
becomes complicated. We thus numerically calculate the Kerr coefficient x'3 = ^(X3) versus 
the pumping Qex and the decoherence rate 75. The result is shown in Fig. [5j We find that 
X3 depends weakly on 75. When 75 becomes non zero, x'z/^ remains to be larger than 1 
MHz. We further calculate the dispersion factor Xi/x'zi the linear absorption factor Xi/x'31 
and the nonlinear absorption factor ^s(x3)/X3 an d plot the results in Figs. [6J [TJ and El We 
see that the unwanted effects increase rapidly with increasing 75. Nevertheless, in the region 
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of current technology 75/27? = 0.5 MHz, these effects are still smaller than the XPM by at 
lease two orders of magnitude. If the decoherence effects can be suppressed to the range 
7 5 /27r = 0.1 MHz in the future experiments, a more pure XPM with negligible dispersion 
and absorption can be achieved. 



V. APPLICATIONS AND DISCUSSION 



As mentioned in the previous section, the molecule induced fluctuation of the TLRs is 
approximately two to three orders lower than the cross Kerr nonlinearity which is on the 
level MHz. In addition, with current technology, the Q factor of a single TLR has already 
been pushed to the order of 10 6 [7|, which yields the decay rate on the level 10 kHz. Since 
both the intrinsic decoherence of the TLRs and the molecule induced decoherence are all 
very small compared with the strength of Kerr nonlinearity, several important inter-TLR 
quantum operations in which only few photons are involved in each TLR could be realized 



with very high fidelities. Recent 
scalable quantum computing [47 



y there are proposals of using the TLRs as qubits to realize 



48| . The TLR states with zero and one photon are used as 



the logical |0) and |1). Therefore, our proposed Kerr nonlinearity can obviously be exploited 
to establish the control phase gates in this system. Another application of the cross Kerr 



nonlinearity is the generation of the inter-TLR macroscopic maximal entangled cat state 



49 



51] . In the first step we tune the molecule decoupled with the TLRs. We then use classical 



pulses to pump the TLRs to the factorized coherent state \a) A \ fi) B and then adiabatically 
tune the molecule to establish XPM with strength X3 between the TLRs. After a time 
t = 7r/x3, the initial \a) A \/3) B evolves to a state 

l*> = \(\<*)a W)b + ha>A W)b + K hP) B ~ \~ a ) A \~P)b) (38) 

Although the coherent states are not orthogonal, the overlap is only 10~ 7 for a mod- 

erate |/3 — a I =4. Superpositions of this kind have no classical counterpart and correspond 
to Schrodinger catlike states. Therefore, the state |^) can be viewed as the macroscopic 
entangled cat state, which can be used as logical elements in quantum computation and can 
offer further insight into the boundary of classical-quantum worlds. We mention that since 
the coherent state spreads over the whole Fock space, such an entangled cat state can hardly 
be prepared by the generalized Eberly-Law algorithm, which has been used to generate the 
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NOON states in the circuit QED system [l7|,ll8|. 



An alternative coupling configuration we can achieve is that the TLR A couples the 
|1) -H- |3) transition while the TLR B couples both the |1) -H- |3) and |2) -H- |4) transitions. 
The coupling between the TLR B and the |1) -H- |3) transition can result from either the 
intention of the experiment or the fabrication errors of the coupling capacitors. In this 
situation, the difference between the mode frequencies plays an important role. When 
ua = wb, an effective Hamiltonian 

Heff = - (39) 

in which the self Kerr effect and the cross Kerr effect co-exist is obtained. When ua 7^ wbi 
since the TLR B is not resonant with the classical pumping field, it can hardly transfer 
population to the state |2), the effective TLR-TLR Hamiltonian thus remains the form 
~9Ai9B2 a ^^ab/(AQ'^ ;x ). From this point of view, it is favorable to use TLRs with different 
eigenfrequencies to realize the XPM, since the difference of their frequencies makes the 
scheme more robust against the fabrication errors. 

Recently theoretical and experimental works have suggested that the large Josephson 
junction can be used to produce SPM and XPM for the TLRs 
a large Josephson junction with the TLRs, one can obtain the Kerr nonlinearity from the 
Taylor expansion of the Josephson energy up to the fourth order. The Kerr nonlinearity 
obtained in this way is smaller than what we have proposed, often on the order of hundreds 
of kHz. Moreover, the critical current noise and the flux noise in the large Josephson 
junction results large fluctuations of the Josephson energy and consequently the severe linear 



4J-|46, 



481 ] . By connecting 



dispersion of the TLRs 
spin-echo like technique 



b eliminate the effect of the fluctuation, one have to use additional 



48] which complicate the quantum gate sequence. Compared with 



these schemes, our alternative XPM scheme may be more robust against the noises and may 
offer more pure Kerr nonlinearity with suppressed frequency drifts. 



VI. CONCLUSION 



In conclusion, in this paper we have shown that an "artificial" multilevel system in cir- 
cuit QED produces the effective XPM between two TLRs with strength much larger than 
previously known. The obtained XPM is very robust against the noises in solid state sys- 
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tern. Compared with the XPM strength, the accompanying dispersion and absorption is 
negligible. Various QIPs can be implemented in this architecture. This work may offer 
improvement to the future scalable quantum computation in superconducting devices. 



Acknowledgments 

We thank Prof. Y. F. Zhu, S. Rebic, Z. W. Zhou, and Y. F. Xiao for fruitful discussions. 
This work is funded by the start funding of HUST, No. 01-24-012018 and No. 01-24- 
012030, the Natural Science Foundation of Hubei province, and the National Natural Science 
Foundation of China, Grant No. 60871018. 



Appendix A: The effective XPM Hamiltonian 

The derivation of Eqs. ( !T0|) and (II ip is provided in this Appendix. As shown in Fig. [Tj, 
taking the decoherence processes |3) — > |1), |3) —> |2), and |4) — > |2) of the atom and the 
decay of the two cavities into account, we write the system-bath Hamiltonian as 

H whole = H sys + H bath , (Al) 

where H damp represents the coupling of the system to reservoir mediating cavity decay and 



spontaneous emission; it takes the form [33 ] 

+oo +oo 

H damp = / ub](u})bj(uj)dw + / up}{w)Pj{w)du> (A2) 

3 =1 > 2 -oo J= 1 ' 2 ' 3 -'oo 
+oo 

+oo 

— oo 
+oo 

— oo 
+oo 



18 



with {bj(u), f3j(u)} the reservoir's annihilation operators at frequency u and {ftj,7j} the 
corresponding decoherence rates. Here we have made the Markovian approximation that 
the system-bath coupling coefficients are constant around the frequencies of interest. 
The master equation of the density matrix p of the system is 

^ = i[p, H sys ] + /tiL[ai]p + n 2 L{a 2 ]p (A3) 
+7iL[ai 3 ]p + 7 2 L[a 23 ]p + 73-^[o r 24]p, 

where L[c\p = 2cpc i — c^cp — pc^c is the Lindbladian form. The expectation value (O) of a 
particular system operator A thus evolves as 

^ = *«>!> = <*>. 

where 

K = i[H sys , 0} + Kl M[ ai }0 + K 2 M[a 2 }0 (A5) 
+ 7l M[a 13 ]0 + 7 2 M[a 23 ]0 + l3 M[a M }0, 

with the M form M[c]0 = 2c^Oc-c^cO -Oc^c. With the help of Eq. (lA"4|l . we calculate the 
time evolution of all the expectation values of the atomic operators. The evolution equations 
of the population operators are 

'" f ' ' - 1 ' 9i (a\ a n + o- 31 a^ + 271 a 33 ^) , (A6) 
g2 \ a\<7M + o"4 2 a 2 j + 27 3 a 4 4^) (A7) 



dt 

d (a 22 ) 
dt 

+ (o"23 + 0-32) + 27 2 cr 33 ) 



d (0-33) 



I g 1 (a\a 13 + cr^a^ + 271033 (A8) 



dt 

- ({tt c ((X 23 + (T 32 ) + 27 2 CT3 3 }) 



d (0-44) 



while the evolution equations of the coherence operators are 



j# 2 (a\(T 2i + a 42 a 2 ) + 27 3 (T4 4 |^) , (A9) 



- = (tt c ai 3 ) + (gie^a-i + fi^ku) , (A10) 
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<^L = (_ ( 7l + 72 + iS ) a 13 - Q c a 12 ) (All) 
+ (9l (o"33 - <Tn)ai) , 

= (~ (73 + cr 14 ) + (5-10-3401 - ^ 2 o"i2a2) , (A12) 

d (cr 32 ) 



((-7i - 72 + i£) 0-32 + ^ c (0-33 - 0-22)) (A13) 

+ ( 92(12034 - 9ia\ a i2 



dt 



= (~ (73 + iA) a 2 4 + fi c o"34) + (52 (o"44 - o- 22 ) a 2 ) , (A14) 
d(au) 

- ((gia\au + 920320-2 



df (-(7i + 72 + 73 -^ + ^A)a- 34 -fi c (T 2 4) (A15) 



When the adiabatic conditions in Eq. (j^J) are fulfilled, the degrees of freedom of the 
atom follow those of the cavities. Therefore, we set the right side of Eqs. (1A6l) - (1A15j) to 
be zero and represent the atomic operators by the annihilation and creation operators of 
the cavities. Suppose initially the atom is prepared in its ground state |1), we can choose 
(°"n) = 1 as a meaningful start and solve the Eqs. (1A6|) - (1A15|) iteratively. The stationary 
values up to the third order of the iteration are 

2 2 
(a n ) s = 1 - ^4°i> (022)5 = ^ a i a i' ( A16 ) 

(0-33)5 = (0u) s = °> 

/ \ 91 01 | , (7l + 72 + 0102 t /Ai7\ 

(<r*)s = -^a, + ^aWa, + ^ + - A) ^3 ( A17 ) 

/ v 0102 f / v 0102 

(73 + iA) 12 2 . (73 + lA) il c 

2 

(032)3 = (0-34)5 = 0, (024)3 = ~ ^ T~^j ^2 a i QlQ2 " 

Setting 73 = we get Eqs. [TU] and dU from Eq. IAT71 

In this Appendix, we have treated a relative simple case. Only few decoherence channels 
are involved, while the atom-photon coupling configuration is also very "clean". For the 
more complicated situations discussed in the manuscript, this systematic method is still 
valid: We first modify H who i e according to the problems we consider and then re-derive the 

20 



evolution equations for the atomic operators; After performing the adiabatical elimination 
and iteration, we can get the stationary values of the atomic operators which contain the 
information of the effective evolution of the two cavity modes. 
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